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A High Signal to Noise Map of the Sunyaev-Zel'dovich 
increment at 1.1 mm wavelength in Abell 1835 
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ABSTRACT 

We present an analysis of an 8 arcminute diameter map of the area around the galaxy 
cluster Abell 1835 from jiggle map observations at a wavelength of 1.1 mm using 
the Bolometric Camera (Bolocam) mounted on the Caltech Submillimeter Observa- 
tory (CSO). The data is well described by a model including an extended Sunyaev- 
Zel'dovich (SZ) signal from the cluster gas plus emission from two bright background 
submm galaxies magnified by the gravitational lensing of the cluster. The best-fit 
values for the central Compton value for the cluster and the fluxes of the two main 
point sources in the field: SMM J140104+0252, and SMM J14009+0252 are found 
to be y = (4.34 ± 0.52 ± 0.69) x 10" 4 , 6.5±2.0 ± 0.7 mJy and 11.3±1.9 ± 1.1 mJy, 
where the first error represents the statistical measurement error and the second error 
represents the estimated systematic error in the result. This measurement assumes 
the presence of dust emission from the cluster's central cD galaxy of 1.8 ± 0.5 mJy, 
based on higher frequency observations of Abell 1835. The cluster image represents 
one of the highest-significance SZ detections of a cluster in the positive region of the 
thermal SZ spectrum to date. The inferred central intensity is compared to other SZ 
measurements of Abell 1835 and this collection of results is used to obtain values for 
y = (3.60 ± 0.24) x 10~ 4 and the cluster peculiar velocity v z = -226 ± 275 km/s. 

Key words: infrared: galaxies - galaxies: clusters: general - methods: data analysis 



1 INTRODUCTION 



The Sunyaev ZePdovich (SZ) effect (|Sunveav fe Zel'dovichl 
1 197(3 ) is the redistribution of energy in the Cosmic Mi- 
crowave Background (CMB) spectrum due to interations 
between CMB photons and hot electrons along the line of 
sight between the surface of last scattering and an observer. 
The main source for the SZ effect is from the hot gas that 
exists in the intra-cluster medium (ICM) of massive galaxy 
clusters llBirkinshawl 19991: ICarlstrom. Holder fe Reesj|2002l : 
iRephaeli. Sadeh fe Shimon! 1 20061 '). The SZ actually com- 
prises two effects. The thermal effect consists of a dimming 
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or decrement in the apparent brightness of the CMB to- 
wards a galaxy cluster at low frequencies and a correspond- 
ing brightening or increment at high frequencies with the 
null crossover point at approximately 215 GHz. The kine- 
matic effect has the spectral dependence of a standard tem- 
perature shift in the CMB which has the same sign at all 
frequencies. The two effects can therefore be distinguished 
from each other with measurements at multiple frequencies. 
Measurements of the amplitude of the SZ thermal 
distortion towards a cluster can be combined with mea- 
surements of the X-ray emission to determine the angu- 
lar diam eter distance, cLa, whose value depe nds on cos- 
mology IjBirkinshaw. Hughes fe Arnaudl Il99lf) . Estimates 
of the Hubble constant have been ma de usi ng this 
technique for a number of c l uster s (e.g. IJonesI 
ICraingdlfl996T ): lHolzapfel et al.l(|l997l ); lTsuboi et"al] 



1995) 



1998 
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iMauskopf etU] J200d); 
(2003); IUdomprasert et al 



Reese! l|2003l ); iBattistelli et~ai1 



(2004J)) and can be used to con 
strain cosmological models. In addition, because the SZ 
surface brightness for a cluster with a given mass is al- 
most independent of redshift, SZ surveys can give infor- 
mation about the evolution of the number counts of clus- 
ters vs. redshift. This depends strongly on the evolution 
of the so-called dark energy or cosmological constant and 
therefore SZ surveys have been identifi ed as one of t he ke y 
probes of the nature of da r k energ y ( e.g. Diego et al.l 1 2002); 
IWeller. Bettve fc Kneissll J2002): iDeDeo. Spergel fc Trac 



2005); lAlbrecht et al.l J2006I ): fehattacharva fc Kosowskv 



2007)). A number of dedicated SZ surveys are already 



producing results, in particular the Atacama Cosmology 
Telescope (ACT) and the South Po l ar Telescope (SP T) 
|Plagge et alj|2009l ; Iffincks et al]|2009l ; iHigh et aljlioiol ). 

Most of the SZ detections reported to date have been 
made at low frequencies corresponding to the SZ decre- 
ment. Follow-up photometric or spectroscopic measurements 
of known clusters at higher frequencies would significantly 
improve the precision in the measurement of the kinematic 
SZ effect as well as constraining possible contamination in 
the low frequency data. Measurements at high frequencies 
corresponding to the SZ increment suffer from confusion 
from emission from dusty galaxies, including background 
high redshift galaxies amplified by the gravitational lensing 
of the cluster as well as increased atmospheric contamina- 
tion from ground-based telescopes. Accurate measurement 
of the SZ increment requires a combination of angular resolu- 
tion sufficient to resolve and remove point sources combined 
with high sensitivity and control of systematics necessary to 
detect the more diffuse SZ signal. 

This paper presents analysis of observations of the 
galaxy cluster Abell 1835 using the Bolometric Camera 
(Bolocam) mounted on the Caltech Submillimeter Observa- 
tory (CSO), situated on the summit of Mauna Kea, Hawaii. 
We also compare the results of this analysis with other sets 
of data taken at different wavelengths. Abell 1835 is one of 
the most luminous clusters observed in the ROSAT cata- 
logue and is well-known as a cooling core cluster. It is also 
known to contain two lensed sub-mm point so urces, SMM 
. 11400 9 +0252 and SMM J14 0104+0252 (see e.g. HvIson et all 
l|2000l ); IZemcovet all l|2007l )). 

The paper is organized as follows: Section 2 describes 
the observations with the Bolocam instrument; Section 3 
describes the analysis pipeline developed for processing the 
data; Section 4 describes the modelling of the data and the 
determination of the characteristic parameters of the model 
(as well as the errors in their values), and Section 5 discusses 
the results of the analysis of the Bolocam data and their 
combination with other literature results. 



2 OBSERVATIONS 

The data used in this work represents approximately 12.5 
hours of observations taken over the course of five nights 
at the end of January/ beginning of February 2006. The 
observations were made with Bolocam operating at 1.1 mm 
(275 GHz). 

Bolocam consists of an array of 105 operational 
neutron-transmutation-doped (NTD) Germanium spider- 



web bolometers, capable of observin g at 1.1 mm or 2.1 mm 
l|Glenn et al.lll998l ; lHaig et al.ll2004) . The system is cooled 
to 270 mK to allow the array to operate close to the photon 
background limit at 1.1 mm. The Bolocam array has a field 
of view of approximately 8 arcmin, and the beam FWHM is 
- 30" at 1.1 mm (and ~ 60" at 2.1 mm). 

The detector is mounted at the Cassegrain focus of the 
Leighton telescope - the 10.4 m CSO dish. The principal 
science targets for the instrument include star forming re- 
gions in the galaxy, blank-field surveys for dusty extragalac- 
tic point sources; blank field SZ cluster surveys, and pointed 
observations of galaxy clusters. 

Observations by Bolocam incorporate a number of dif- 
ferent scan strategies. The majority of observations use 
raster scanning or lissajous scanning where the entire tele- 
scope is constantly in motion modulating the array position 
on the sky. 

The data presented in this paper was taken using jiggle 
mapping. Jiggle mapping uses a combination of 'chopping' 
the secondary mirror during observations, i.e. switching the 
secondary from side-to-side, and 'nodding' the telescope, i.e. 
changing the position of the telescope during a scan. The 
beam moves from being 'on-source' (imaging the region of 
the target and chopping away from it), to 'off-source' (imag- 
ing a nearby region of blank sky and chopping onto the tar- 
get), then returns to being on-source once again (see Fig. 1). 

Jiggle mapping is a means of removing sky background 
from data, and works by differencing signal from so-called 
'target' and 'reference' beams (centered on the science target 
and a region of 'blank' sky, respectively). These beams are 
also referred to as 'on' and 'off' beams. While jiggle-mapping 
is an efficient method of subtracting sky-noise, care has to 
be taken to ensure that the chop 'throw' (i.e. the maximum 
angular dispacement of the secondary during the chop) is 
large enough that the off beam does not pick up signal from 
the source itself (subject to mechanical limitations). For this 
reason, jiggle-mapping is not always suitable for imaging 
extended sources, but is good for observing point sources. 

The chop frequency for this data set was 2.25 Hz, giving 
a chop period of ~ 0.44 s (compared to a sampling rate of 
0.02 s). The chop throw was 90", while the displacement 
between on and off beams during the nods was 2.5'. The 
telescope typically remained in the 'on' or 'off' phase for 
periods of approximately 10 seconds. 

The full set of data consisted of fifteen individual sets 
of observations of approximately 50 minutes each. The ob- 
servations were separated by shorter observations of planets 
and other sources with known position and flux, which could 
be used for the calibration and pointing. 



3 PIPELINE DEVELOPMENT 

Because jiggle-mapping is a relatively new strategy for Bolo- 
cam, there was initially no pipeline in place to analyse the 
raw data. An independent pipeline was written specifically 
for this observing mode. 

In its original form, the time-stream data had been 
sliced (i.e. separated into individual data-vectors represent- 
ing the signal from individual bolometers during specific ob- 
servations), but no other processing had taken place. The 
pipeline for jiggle data includes the following elements: 
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Figure 1. Nodding and acquisition pattern (top). A nod position 
of +1 implies the telescope is on-source, whereas nod position -1 
is when the telescope is off-source. The pattern is a combination 
of acquisition and positional data, so that where the pattern falls 
to zero this reflects where the telescope is not acquiring data 
(usually while the dish is moving to a different location on the 
sky). A sample of the raw data demonstrating the chopping is 
also given (bottom). 



1. Initial cleaning to remove obvious sources of noise 
using average subtraction; 

2. Identifying the nodding sequence; 

3. Identifying phase differences between the signal from 
the chop control and the data (due, for example, to differ- 
ences in timing between the telescope clock and those in the 
telescope control terminals); 

4. Deconvolving the data for each nod position ('remov- 
ing' the chopping); 

5. Extracting the signal (differencing the deconvolved 
signal from the nod positions); 

6. Calculating errors for the signals; 

7. Determining pointing corrections for each map and 
hence reconstructing the directional information for each 
map; 

8. Calibrating the maps; 

9. Saving the deconvolved datasets, 

10. Coadding maps from individual observations. 



3.1 Deconvolving 

We deconvolve the chopped data using a fit to the time 
stream signal monitoring the chopper position corrected for 
a phase difference between the chopper data and bolometer 
data. If the signal being detected by an individual bolometer 
is modulated by the chopping according to: 



F(i) = F s sin[v c t + <p c ] 



(1) 



where P(t) is the time(i)-varying flux received by the 
bolometer, F s is the flux of the source, i/ c is the chop fre- 
quency, <f) c is the chop phase, one can then multiply this by 
a fit and integrate according to: 



S(t) 



/F(t)sin[i/ft + f ]dt 
f sm 2 [vft + <f>f]dt 



(2) 



where i>{ is the fit frequency, <j>{ is the fit phase. 

With an accurate enough fit (ut ~ i/ c , <f>{ ~ </> c ), there- 
fore, S(t) —¥ F s . The value of V{ is measured from a fit to 
the chopper data, and the integration needs to be over a 
complete number of chopping cycles to avoid spurious sig- 
nal being introduced into the maps. 



3.2 Pointing 

The individual detector pointing data was reconstructed 
in three stages. First a map was made of the deconvolved 
data. The directional information was taken from the time- 
streams of azimuth and elevation stored with the raw data 
(converted to ra and dec) and previously determined offsets 
(to convert the values stored at the telescope with the actual 
azimuth and elevation of the centre of the telescope beam) 
were added. The mean of these values were compared to the 
recorded position of the source to obtain a 'mean' offset. 

The directional information for the map was then re- 
constructed to include the mean offsets. The location of the 
centre of emission for the source in the map was determined 
using a simple \ 2 n ^ °f the data immediately around the 
cluster to a gaussian. This position was then compared once 
more with the source values to determine a 'fine' offset. 

The directional information was then undated again to 
include both the mean and fine offsets. 



3.3 Calibration 

The calibration for Bolocam is determined using the resis- 
tances of the bolometers as determined by the dc level of 
the lock-in amplifier output signal, which is expected to be 
directly rela ted to the flux calib ration for a set of obser- 
vations (see lLaurent et al.l (|2005f )). If observations can be 
made of a number of sources of known flux at different lev- 
els of sky loading a plot can be made of the dc level against 
calibration factor. We fit a second-order polynomial to the 
measured calibration points to extend the calibration to the 
full range of observation conditions. In this way, we can de- 
termine the correct calibration for an observation within an 
appropriate range of dc level. 

These calibration curves are not expected to change sig- 
nificantly over time. However, given that the previous set 
of data was taken in May 2004, using a different scanning 
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Figure 2. Calibration data. The May 2004 calibration curve is 
shown as the solid line. The error bars in this figure are statistical 
errors from the measured data. They do not, therefore, include 
intrinsic uncertainty in the source brightness. The data does not 
deviate substantially from the May 2004 curve. Legend: Verti- 
cal crosses: 0420m014; open triangles: 0923p392_4cp39.25; open 
diamonds: 1334-127, open squares: 3c371. 



technique, we performed a full set of calibration observa- 
tions roughly every 20 minutes during the science obser- 
vations. These calibration sources were chosen to be close 
to the science observation target and included: 0420m014; 
0923p392_4cp39.25; 1334-127, and 3c371. The resulting cal- 
ibration curve is given in Fig. 2. 

The observations of the calibrators agree reasonably 
well with the May 2004 calibration curve. The point sources 
used in the plot were all secondary sources and as a result 
it was not clear how precise or stable their fluxes would be. 
Given this, the deviations between the old calibration curve 
and the sources were not considered to be great enough to 
warrant revising the calibration values and adopting a new 
functional form for the calibration, and the May 2004 cali- 
bration was used throughout. 



4 MODELLING AND PARAMETER 
ESTIMATION 

Once maps had been made of the individual observations, 
they were coadded to produce a single image of Abell 1835 
which was then also convolved with a Gaussian PSF with 
FWHM ~ 30.6 " , corresponding to the best-fitting Gaussian 
to the Bolocam beam (Fig. 3). 

A model was produced to simulate the observations and 
derive estimates of the fluxes of the sources. Fig. 3 clearly 
shows the two point sources arranged on either side of the 
emission from the cluster (which is assumed to be entirely 
signal from the SZ). These point sources have been de- 
tected previously at 850 /im using the Submillimeter Com- 
mon User Bolometer Array (SCUBA) on the Ja mes Clerk 
M axwell Telescope (JC MT) (jzemcov et al.ll2007T ) as well as 
by llvison et al.l (|2000r). The chara c teristi cs of these point 
sources, reported by IZemcov et al.l (|2007f ) are given in Ta- 
ble 1. 

The SZ emission from the galaxy cluster itself was mod- 



elled assuming that the signal was dominated by the thermal 
SZ and assuming an isothermal beta radial profile: 



r(o p ) = A (i + %. 



(MP) 



(3) 



where T(0 P ) is the flux profile of the cluster as a function 
of projected angle 9 P , 9 C and /3 are fitting parameters A is a 
scaling parameter, defined by: 

.,26 



A = Affloyog(x) x 10 



(4) 



Here: Af2 is the solid angle of the pixels in the (model) map 
(the map was later convolved with the Bolocam beam such 
that the final flux measurements were given in flux/beam), 



(" 



bTqml 



h 



is the blackbody emission of the CMB 
at the redshift of the cluster, yo is the so-called 'Compton 
parameter' at the centre of the cluster (the Compton pa- 
rameter in general, y, can be characterised as: 



y : 



ksT e ,. 
n e a T ir-dl 



(5) 



where n e and T e are, respectively, the electron density and 
electron temperature along the line of sight), and g(x) is 
a function that depends on the dimensionless frequency at 
which the observations are being made, x = 



hi/ 



ksTCMB ' 



sW = x* 



expx 



(expx — 1) 



r (xcoth(a;/2) -4) 



(6) 



The temperature reported in the lite rature for Abell 
1835 varie s between app ro x. 8- 10 keV (e.g. [A llen & Fabian 
(1998); IPeterson et alJ IJ200ll ): iKatavama fc Havashidal 
(2004)). At these temperatures, relativistic corrections to 
the SZ can be on the order of ~ 10% and have to be taken 
into account during the analysis. We u s ed th e corrections 
published by lltoh. Kohvama fc Nozawal (|l998T ) , up to fifth 
order in 8, where: 



k B T e 
m e c 2 



(7) 



An ideal sky model was produced including extended 
emission from the cluster and emission from the point 
sources at positions given bv lZemcov et al.l (|2007T ). The free 
parameters of the model were the cluster geometry (9 C and 
/3); the comptonization (yo), and the fluxes of the two point 
sources J"i and Ti (corresponding to SMM J0140104+0252 
and SMM J14009+0252, respectively). For each set of pa- 
rameters, we produced a high resolution map with pixel size 
of 1" . This map was then convolved with the Bolocam beam 
to produce a 'real-sky' map of the cluster at the resolution 
of the telescope. This in turn was then 'observed' by using 
the data ra and dec values to reproduce the effect of the 
jiggling, thereby producing a simulation of what we expect 
to observe given the assumed values of flux (for the point 
sources) and yo. Finally, the chopped map was rebinned to 
the same pixel size as the raw data maps. By comparing 
the maps corresponding to different values of the parame- 
ters with the data, we construct a multi-parameter likeli- 
hood space, which could be searched to extract best fitting 
models. 

Values for these parameters were determined by search- 
ing for the set that minimized the x 2 value between the 
model map and the data map. 

Generally, SZ observations are not able to constrain 



A High Signal to Noise Ratio Map of the Sunyaev-Zel'dovich increment at 1.1 mm wavelength in Abell 1835 5 






Figure 3. (Top left) S/N map of the raw data, convolved with the Bolocam beam; (top right) 3-parameter fit (/3 = 0.69, 9 C = 33.6") 
S/N map convolved with the Bolocam beam, (bottom) Map of residuals, formed by taking the difference between the convolved model 
and data S/N maps. 



Source id 


ra (hh:mm:ss) 


dec (dd:mm:ss) 


850 fim Flux (mjy) 


450 /an Flux (mjy) 


SMM J14011+0252 
SMM J14009+0252 


14:01:04.7 
14:00:57.5 


+02:52:25 
+02:52:49 


14.6+1.8 
15.6+1.9 


41.9+6.9 
32.7+8.9 



Table 1. Point source positions and fluxes reproduced from llvison et al.l (2000) 



the parameters 8 C and /?, since the resolution of such ex- 
periments is not accurate enough to constrain their val- 
ues to a better degree of accuracy than equivalent X- 
ray (or other wavelength) measurements. Recently, how- 
ever, computer simulations have suggested that using X- 
ray derived profiles for fitting SZ observations can lea d 
to bias in the derived values of y (|Hallman et al.l 120071 ). 
Furthermore, the values of the profile parameters quoted 
in the literatur e can va r y con si derably (see, e . g. va l- 
ues reported in |jia et all J2004); iPeterson et all 1 200 it) : 



1 J200l): 
(|200ll ) : |zi 



ISchmidt, Allen fc Fabianl (|200ll ); izemcov et alt (|2007h ) 



In this work, therefore, a range of approaches have been 
adopted. Fits to y, Ti and T2 were made using two differ- 
ent sets of values for 6 C and /3 obtained using Chandra X-ray 
data and interferometric SZ data from the Owens Valley Ra- 



dio Observatory (OVRO) and th e Berkeley Illino i s Mar yland 
Association (BIMA), taken fror uLaRoque et alt (|2006l ). 

The cores of galaxy clusters can exhibit sharp peaks 
in emission due to the presence of cooling flows, and for 
this reason a single isothermal beta model fit to all the 
X-ray data in a set o f observations is not always suitable. 
lLaRoque et alt (|2006l ) use two different strategies to deal 
with this. In the first the X-ray data set is restricted to ex- 
clude the region within 100 kpc of the X-ray centre, and the 
profile parameters are then obtained by a joint fit to this 
X-ray data and the entire set of SZ data. It is possible to 
use the entire SZ dataset since SZ emission has a weaker 
dependence upon electron density than X-ray observations, 
and the SZ profile is not expected to be significantly affected 
by the presence of cooling flows. (Indeed in general, the SZ 
profile is less sensitive to the detailed physics in the cores of 
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galaxy clusters (|Motl et alj|2005h .) In the second approach 
they employ a double-beta model which has been used to fit 
all the X-ray data in ob servations that include emi ssion from 
cooling flows ( s ee e.g . iMohr. Mathiesen fc Evrardl (1999)). 
lLaRoque et al.l (|2006l ) also fit to the SZ data only (although 
this fit is really only to 9 C since the value of /3 is fixed to the 
value obtained for the isothermal beta X-ray/ SZ fit). 

For the purpose of this analysis the results from the 
isothermal beta fit to the restricted X-ray data and SZ data 
(hereafter LaRoque X-ray+SZ) and the fit to the SZ data 
(hereafter LaRoque SZ only) were used. The values of 9 and 
P they obtain for these fits are, for the combined X-ray and 
SZ data: /? = 0.69, 9 C = 33.6", and for the SZ-only data: 
/3 = 0.70, 9 C = 50.1". 

It was decided, therefore, to fit the data from the Bolo- 
cam observations with five combinations of profile parame- 
ters: = 0.69, 9 C = 33.6"; /3 = 0.69, 9 C = 50.1"; /3 = 0.69; 
9 = 33.6", and 9 = 50.1", where the other parameters are 
fitted for in each situation. 

The goodness of fit for a particular set of parameters 
was evaluated using a \ 2 test. The value for \ 2 f° r a given 
model was defined simply as: 



2 

X = 



E 



(datai — modeli) 



(8) 



where data^ denotes the signal in the i'th pixel of the 
data map, modeh represents the signal in the i'th pixel of 
the model map, and Oi is the noise in the i'th pixel of the 
noise map. The best fitting parameters for a given analysis 
run were, therefore, those that minimized this value of \ 2 '. 

For the cases in which there were four free parameters 
to be searched through, the individual searches took a long 
time to process. The general method of locating a minimum 
was to sample the parameter space between 'reasonable' lim- 
its, locate a minimum and examine variation in the \ 2 val- 
ues for each parameter while keeping the other parameters 
fixed. This variation was in general smooth and reasonably 
quadratic, so that the first derivatives of the variation was 
typically linear. This could easily be fit and an estimate 
made of where the first derivative was zero. The limits of the 
search were then reset to centre on the new approximation 
for the minimum and the process repeated until the new es- 
timates of the minima were similar (within some tolerance) 
to the previous ones. This approach allowed the minimum 
values to be located within a reasonable amount of time. 

Errors on the values of the fit parameters were found by 
deriving the inverse of the Fisher matrix for the observations 
and evaluating its diagonal elements. For parameters with a 
gaussian distribution, the Fisher matrix Fij is given by: 



1 



2 \ dpidpj 



(9) 



where p t represents the vector of parameters being fit, and 
the second-order differential is evaluated at the value of p at 
which the \ 2 value is minimized. 

This kind of analysis relies on there being no covari- 
ance between different pixels in the data maps. The level of 
covariance is expected to be low due to the chopping, and 
because the only cleaning of the data involves an average 
subtraction. This should introduce less pixel to pixel covari- 
ance than other cleaning methods. Jacknife maps were also 



made to check these assumptions, and were found to show 
no evidence for significant levels of covariance. 

The components of Fij were evaluated using a finite dif- 
ference scheme and arrays of \ 2 values around the minimum 
values. The \ 2 values were obtained by varying the fit pa- 
rameters as well as the position of the central emission. The 
results of this analysis are shown in Table 2. 



5 RESULTS 

The results of Table 2 shows that the 4-parameter fit with 
9 C — 50.1" gives the best (\ 2 ) value. The best fit values of 
the central Compton parameter as well as the fluxes of the 
point sources are, therefore, found to be: yo = (4.19±0.48) x 
10" 4 ; 7.17±1.84 mjy, and 12.18±1.77 mjy respectively. The 
errors on the values of j/o in Table 2 include statistical er- 
rors as well as errors due to pointing and uncertainty in the 
model parameters. They do not include systematic errors 
(e.g. due to dust contamination, kinetic SZ, etc.). The sys- 
tematic uncertainties are discussed in greater detail below. 
The value of /? for our best fit model is substantially 
different from values quoted by other groups, who use X- 
ray data with much better resolution to calculate their pro- 
file parameters. It is common practice to treat these mea- 
surements as more reli able in determining cluster profile pa- 
rameters (although see lHallman et al.l l)2007ri ). Furthermore, 
since a comparison of our data to that of other groups was 
planned, the cluster model for which intensity measurements 
were being compared had to be the same. For these reasons, 
it was decided that, in forming the comparison with data at 
other frequencies, the results for the more 'canonical' values 
the profile parameters should be used. Of the two fits with 
more typical profile parameters, the /3 — 0.69, 9 C — 33.6" 
fit has the lower \ 2 value, and is the model used for the 
spectral fitting. 



ADDITIONAL SOURCES IN THE FIELD 

In order to reliably simulate a series of chopped observa- 
tions, all sources in the field need to be accounted for. The 
j3 = 0.69, 9 C — 33.6" fit model is also shown in Fig. 3, 
convolved with the Bolocam beam. The final map in Fig. 3 
displays the difference between the model map and the raw 
data, convolved with the Bolocam beam to emphasize any 
differences between them. The difference map appears to 
indicate that the cluster signal and the point sources have 
been removed effectively and that, therefore, the model is a 
reasonable simulation of the data. The signal at the edge of 
the field is mostly spurious. It could be suggested that there 
is a ~ 3-ct detection (in the raw map) of a point source just 
below centre of the image field, but there are no obvious 
listed sources located at the position of the candidate detec- 
tion, and it appears that most likely to be a spurious feature 
of the noise. 

Using noise maps calculated in the analysis, it was pos- 
sible to form a histogram of the difference between the S/N 
data map and the simulated map (see Fig. 4, where the 
data is plotted on logarithmic axes to determine whether 
there was significant deviation between the residuals and 
the fit toward the limits of the data). The Gaussian form 
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Table 2. Parameter estimates (1 mm Bolocam observations) and x 2 values; Number of degrees of freedom 
= 3013. (p) denotes that the model included a 1.8 mjy point source at the cluster centre. 




0.003 F 



Residual S/N 

Figure 4. Signal to noise histogram of the residual map obtained 
by subtracting the best fit model from the data, plotted on loga- 
rithmic scales. The Gaussian best fit to the histogram (the solid 
line in the log plot) has cr = 1.01 



of the histogram (quadratic in Fig. 4) is consistent with the 
residual signal being dominated by noise, rather than other 
point sources. The histogram was fit with a Gaussian, and 
the best fit has a standard deviation close to 1.0, indicating 
that the noise had been correctly estimated. There may be 
some suggestion of deviation in the high S/N region between 
the histogram and the fit, but this is attributed to regions 
of lower coverage in the outer sections of the map, and not 
believed to represent a real excess of sources. 



DUST CONTRIBUTION 

At longer wavelengths the contribution from dust emission 
is expected to become more important and is a source of 
confusion. The literature does not provide an entirely con- 
sistent view of whether or not Abell 1835 has significant dust 
emission. While the cluster's brightest cluster galaxy (BCG) 
exhibits significant CO emission which is believed to trace 
dust, it would be expected that the dust emission would be 
extremely bright at 450 /im. There are a limited number of 
results in th e literature that re port emi ssion from Abell 1835 
at 450 fim. Ilvison et al.1 l|2000l ) and lEdge et all (QUI) re- 
port strong emission, consistent with a signifi cant dust con- 
tribution, whereas the 450 /im result given in lZemcov et al.l 



0.001 r 



-0.001 



-0.002 b 




Figure 5. The SZ spectrum of Abell 1835. The 30 GHz re- 
sults (open triangle, x ~ 0.5) are taken from.LaRoque, et al. 
(2006). The points at 142 and 221 GHz (open diamonds) are 
from Mauskopf, et al. (2000), while the open square represents 
the results reported here. All points have been renormalised to 
the central intensity of the SZ effect for a cluster model with 
8 = 33.6" and = 0.69. The line is the best fit SZ spectrum 
with yo = 3.60 X 10~ 4 and v z = -226 km/s. 



(2007) of -2 ±13 mjy is consistent with there being no emis- 
sion from the cluster at this wavelength. 

Assuming a 'wor st case' scenario , h owever, in which the 
emission reported by Ilvison et al.l (J2000J) at 450 and 850 fim 
are assumed to be correct, and all the signal at these wave- 
lengths is attributable to the cD galaxy, a model can be 
formed of the dust emission using a greybody spectrum, as- 
suming a dust temperature ~ 30 K and an emissivity which 
varies as v . The flux of the cD galaxy at 1.1mm using this 
model is expected to be ~ 1.8±0.5 mjy. This value is consis- 
te nt with the SEP pro duced for Abell 1835's BCG reported 
bv lEgami et al.l (|2006f ). Introducing a point source into the 
cluster model at the X-ray centre and then marginalizing 
again over values of yo, we obtain a best fit value for yo of 
(4.34 ± 0.52) x 10" 4 , with a X 2 value of ~ 3062. 



OTHER SOURCES OF ERROR 

Along with the dust contribution, it was important to ac- 
count for other sources of error in the result. These are dis- 
cussed individually below: 
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Systematic errors 



Error (as % of result) 



Calibration 
Kinematic effect 
CMB confusion 
Dusty galaxies 
Dust emission 
Effective temperature 
Uncertainties in 
King model parameters 
Others (clumping, etc) 



11.7 
9.0 
1.0 
5.4 
7.2 

+1.5 
-2.6 

0.6 



Pointing and model 

error (incl. point 

source 

Total systematic 

Total error 



0.52 mjy 



0.69 mjy 
0.87 mjy 



Table 3. Error budget for Abell 1835 observations 



Calibration 

We estimated the calibration error from the noise-weighted 
dispersion of the measured point fluxes relative to the May, 
2004 model. This gives a calibration error of ~ 10.6%, which, 
when combined with a 5% error due to uncertainties in the 
Mars model that all the Bolocam fluxes are dependent on 
gave an overall calibration error of ~ 11.7%. 



Kinematic effect 

The kinematic SZ effect (see below) is due to the bulk motion 
of a cluster. In general, the kinematic effect is much less 
significant than the thermal effect that has been discussed so 
far. The bulk motion of clusters in a concordance cosmology 
predicts an rms pecul iar velocity of around 300 km/s, but 
iMauskopf et al.l (J2000j) obtain an estimate of the velocity of 
Abell 1835 of around 500 km/s. A velocity this high at 273 
GHz produces a kinematic effect that is approxiamtely 9% 
of the thermal effect. 



Confusion sources 

The main sources of confusion in observations of the SZ ef- 
fect include the CMB background and dusty galaxies. Typ- 
ical CMB temperature anisotropy signals based on the cur- 
rent concordance model spectrum was simulated and 'ob- 
served' using our scan strategy. These produced an rms 
equivalent to a 1.5% error in the measured cluster signal. 
Dusty galaxies at 1.1 m m contribute an rms signal of ap- 
proximately 0.5 mjy (e.g. iBlainl 1)19981 )). which corresponds 
to a 5.4% error in the measured cluster signal. 



Physical model uncertainties 

Other sources of error related to the parameters of the phys- 
ical model of Abell 1835 itself. Uncertainties in the electron 
gas temperature of 10 - 15 % cor respond to an e r ror in the 
estimate of y a of up to ~ 2.5%. iLaRocme et al.l (|2006f l re- 
ports errors on the model parameters d c and /3 of ±1.0" n 
and ~ 0.01, respectively. When these variations are intro- 
duced into the model to see what effect they have on the 



value of j/o obtained by the fit, their effect is found to be 
small, and represent an uncertainty on the level of 0.6 % 
on the final result. Other potential sources of error include 
(e.g.) clumping in the cluster gas, but the SZ effect is rela- 
tively insensitive to the detailed physics of cluster gas, these 
effects were estimated as contributing no more than ~ 2% 
to the final result. 

A breakdown of the error budget for the observations 
presented here is given in Table 3. 



To ESTIMATES 

The value of j/o given by the 4 parameter fit means that the 
observations detect the cluster with a significance of 10.5a. 
This is one of the highest significance detections of a clus- 
ter in the positive region of the SZ spectrum to date. A 
plot of the SZ spectrum of Abell 1835 based upon these val- 
ues is given in Fig. 5. These values have were fitted with 
a spectrum that included kinetic and thermal effects (non- 
relativistic as well as relativistic). The kinetic effect follows 
the form: 



^kinetic — TePll\X) 



(10) 



where r is the optical depth of the cluster gas, given by 
j/o I /" — ) I , /3 = ^ (where v z is the peculiar velocity of the 

galaxy cluster), and h(x) = x 4 -. — cxp f..,y ■ The kinetic effect 
dominates the SZ spectrum around the null point (v ~ 217 
GHz, x ~ 1.9), but makes only a small contribution to the 
spectrum at frequencies above or below this. 

The free parameters of the fit to the SZ spectrum were, 
therefore, j/o (which characterizes the thermal SZ), and v z 
(which characterizes the kinetic effect). The values of these 
parameters that optimized the fit were found to be: j/o = 
(3.60 ± 0.24) x 10~ 4 , and v z = -226 ± 275 km/s. This is 
consistent with prev i ous measuremen t s of v z for Abell 1835 
iBenson et all i|2004f ); IMauskopf et all (|2000r i and represents 
the most sensitive test for peculiar velocity in an individual 
galaxy cluster using the SZ effect to date. 
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6 CONCLUSIONS 

The observation and analysis of 1.1 mm Bolocam jiggle- 
map data of Abell 1835, including details of the process 
of converting the raw data to maps of the cluster, have 
been discussed. A parameter search has been used to deter- 
mine the best fit values of the central Compton parameter 
for the cluster and the fluxes for the point sources SMM 
J140104+0252 and SMM J14009+0252. These values are 
found to be y = (4.68 ± 0.48 ± 0.82) x 10" 4 , 6.5±2.0 mjy 
and 11.3±1.9 mjy, respectively with a statistical signal-to- 
noise of 10.5, 3.3 and 6.0 respectively. If the model is refined 
further to include a point source, representing dust emission 
from the cD galaxy in Abell 1835, the Compton parameter 
for the cluster is found to be y = (4.34 ±0.52 ±0.69) x 10 -4 
The value for yo was compared to other literature results 
and found to agree well. The SZ spectrum for Abell 1835, 
based upon these values for j/o was fit with the full SZ form, 
including both thermal and kinetic effects and relativistic 
corrections up to 7th order. The values of yo and v z which 
optimized this fit were found to be (3.60 ± 0.24) x 10 -4 and 
—226 ± 275 km/s, which are both in agreement with previ- 
ous literature results. It is also concluded that, in order to 
better evaluate the contamination from dust and accurately 
characterize the spectrum of point sources in the same field 
as SZ clusters more results need to be obtained at shorter 
wavelengths in order to obtain more a detailed spectral en- 
ergy distribution. 
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